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The recent measurement of negative higher-order Kerr effect (HOKE) terms in gases has given rise to a 
controversial debate, fed by its impact on short laser pulse propagation. By comparing the experimentally 
measured yield of the third and fifth harmonics, with both an analytical and a full comprehensive numerical 
propagation model, we confirm the absolute and relative values of the reported HOKE indices. © 2011 
Optical Society of America 
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In a recent experiment, we have shown that the elec- 
tronic optical Kerr effect in Ar, N2, O2, and air exhibits a 
highly nonlinear behavior versus the applied intensity [1] , 
resulting in a saturation of the nonlinear refractive in- 
dex observed at moderate intensity, followed by a sign 
inversion at higher laser intensity. This observation has 
a substantial impact on the propagation of ultrashort 
and ultra-intense laser pulses, especially in the context of 
laser filamentation [2-5] where the higher-order Kerr ef- 
fect (HOKE), rather than the defocusing contribution of 
the free electrons, can play a key role in the self- guiding 
process [6], especially at long wavelengths [7] and for 
short pulses [8]. However, this issue is still controver- 
sial [9-11]. Therefore, an independent confirmation of 
our measurement of the HOKE is still needed. Recently, 
Kolesik et al. [9] have proposed such test, based on the 
comparison of the yields of the third harmonic (TH) 
and the fifth harmonic (FH) radiations generated by the 
nonlinear frequency up-convcrsion of a short and intense 
laser pulse in air. Based on numerical simulations, they 
suggested that, considering the HOKE indices, "the rela- 
tive strength of the FH to the TH should reach values 
of the order of 10 _1 " while, if omitting them, "this ratio 
should be about 4-5 orders smaller" [9]. 

So far, no measurement of the yield of FH versus the 
TH have been achieved in air. However, Kosma et al. 
[12] measured the yields of TH and FH produced by 
a short and intense laser pulse in argon. The present 
paper aims at confronting the results of this experiment 
to predictions based on the HOKE in argon [1]. 

In the first part, we confirm the ratio of the recently 
measured non-linear indices [1] based on the analytical 
description of the harmonic generation. In the second 
part, a comprehensive model including linear and non- 
linear propagation effects such as dispersion, self-phase 
modulation, ionization, and Kerr effect, is presented. 

For a focused laser beam propagating linearly, the har- 
monic power of the qth harmonic in the perturbative 



regime is given by 

P q = A q N 2 \J q (bAk)\ 2 , 

where TV is the atomic density of the medium and 
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with Pi, ui, and wo the power, the angular frequency, 
and the beam waist of the incident beam, respectively 
[13,14]. x is t ne gth-order microscopic nonlinear sus- 
ceptibility (q = 3, 5) given in SI units, n,j are the linear 
refractive indices at the fundamental (j = 1) and har- 
monic frequencies (j = 3,5), eo is the permittivity of 
vacuum, and c is the speed of light. J q is a dimcnsionless 
function that accounts for the phase matching 
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match, with n q — n\ proportional to the pressure, and 
kj (J = 1, q) the wave vectors, b the confocal parameter, 
L the length of the static cell, and / the position of the 
focus with respect to the entrance of the static cell [15]. 
According to Eqs. (1) and (2), the ratio of the FH to the 
TH power is 
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where have been approximated to unity in Eq. 2. A3 
and A5 refer to the different atomic densities at the pres- 
sures maximizing the harmonic conversion for the 3rd 
and 5th orders, respectively. This equation provides a 
direct relationship between the power ratio of the har- 
monics and the ratio of the corresponding non-linear sus- 
ceptibilities. The latter are related to the nonlinear re- 
fractive indices through the relation [7] 
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Fig. 1. (Color online) Analytical calculation of the pres- 
sure dependence of the 3rd (solid blue line) and 5th 
(dashed red line) harmonics in argon. 
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In the experiment by Kosma el al., 6=7.8 cm, 
wo = 100 fim, L=1.8cm, f=L/2, and the wavelength 
Ai=810nm [12]. The fundamental power, calculated 
from the input energy Ei=710 p3 and the pulse du- 
ration Ti=12fs, is Pi = 59 GW. They observed that 
the pressure maximizing the TH power ranged between 
160 mbar [12] and 250 mbar [16], for similar experimen- 
tal conditions. One single maximum, around 50 mbar, 
is observed for the FH. The maximum energies of the 
TH and FH measured at the respective optimal pres- 
sures reported in [12] arc 140 and 4nJ, respectively, 
while the pulse duration was estimated to be llfs for 
both harmonics [12, 16]. This leads to a power ratio 
P5/P3 = 0.028. According to Eq. (6), where J q of Eq. (3) 
has been calculated using n\ = 1.00028, n| = 1.00030, 
and rig = 1.00035 for the values of the refractive in- 
dex of argon at 1 bar at 810, 270, and 162 nm, respec- 
tively [17], the corresponding ratio of the HOKE indices 
is |n4/n2| = 6.8 x 10~ 19 m 2 /W. This value confirms, 
within a factor of 2 compatible with the experimen- 
tal error, the ratio of the experimental HOKE indices 
7i 2 = lO^W/W 1 and 774 = -3.6 x 10~ 42 m 4 /W 2 [1], 
resulting in I714/712I = 3.6 x 10 _19 m 2 /W. The agreement 
is remarkable, especially considering the simplicity of the 
analytical model used. 

Further comparison with the experiment was per- 
formed by computing the value of iV 2 | J q \ 2 as a function 
of the argon pressure relying on Eq. (3) (Fig. 1). This 
function should reflect the pressure dependence of the 
harmonic powers. The analytical model predicts a max- 
imum at about 300 mbar for TH, in line with the exper- 
imental results. It yields three maxima between and 
400 mbar for FH, the first of them close to the observed 
optimum pressure for the FH. This oscillatory structure, 
which is due to the periodic phase matching, was not 
observed in the experiments [18] probably due to non- 
linear propagation effects which are not considered in the 
analytical model. 



To overcome these limitations and take into account 
the perturbations of the fundamental pulse during its 
propagation through the gas sample, as well as the ef- 
fect of the HOKE indices on the phase matching, we have 
solved the unidirectional pulse propagation equation for 
the experimental conditions of Kosma et al. More pre- 
cisely, assuming a cylindrical symmetry around the prop- 
agation axis z, the angularly resolved spectrum E(kj_, lo) 
of the real electric field E(r, t) follows the equation [19] 

where k z = ^k 2 {u) — Pnl (resp., J) is the angu- 
larly resolved nonlinear polarization (rcsp., free charge 
induced current) spectrum, and k{u) = n ^' UJ . 

The nonlinear polarization Pnl is evaluated in the 
time domain as P NL = x (3) £ 3 + x (5) £ 5 + x {7) E 7 + 
X^E 9 + x^^E 11 . Since the nonlinear polarization is 
defined from the real electric field, Eq. 7 captures with- 
out any modifications all frequency-mixing processes in- 
duced by the total field. For numerical stability concerns, 
we considered only the part responsible for the refractive 
index change around ujq, neglecting harmonics genera- 
tion induced by the terms proportional to E" 7 , E 9 , and 
E 11 . The current induced by the free charges is calcu- 
lated in the frequency domain as J = ^~ ^X+J 1 P £ ' wnere 
e (resp., m e ) is the electron charge (resp., mass), v e is 
the effective collisional frequency, and p is the electron 
density which is evaluated as 

d tP = W(I){ Pat - P ) + ?-I~p P \ (8) 

where W(I) is the ionization probability evaluated with 
the Keldysh-PPT (Perelomov, Popov, Terent'ev) model 
[3], pat is the atomic number density, a is the inverse 
Bremsstrahlung cross-section, /3 is the recombination 
constant (negligible on the time scale investigated in the 
present work), and / is proportional to the time- averaged 
<£ 2 ). 

Figure 2 displays the harmonics intensity as a func- 
tion of argon pressure for an input pulse and a detec- 
tion geometry matching the experimental parameters: 12 
fs pulse duration (FWHM) , 700 p J input energy, and a 
beam radius of 4 mm before focusing. In order to mimic 
the experiment, the pulse first propagates in vacuum up 
to the position of the cell (99.1 cm after the /=1 m lens). 
After this focusing step, the pulse propagates over 1.8 cm 
in the argon cell. The optimal pressure for the FH is 50 
mbar, in full agreement with the experiment [12]. The 
reduction of the second and third maxima of the FH, as 
compared to Fig. 1, results from the phase mismatch in- 
troduced by the HOKE at large pressure. The TH yield 
is maximal at 260 mbar, similar to the value reported 
in [16]. In full agreement with the experiment by Kosma 
al. [12], the ratio at 50 mbar is about 0.1 and becomes 
even larger at reduced pressures. Furthermore, the total 
FH and TH energies at their respective optimum pres- 
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Fig. 2. (Color online) Numerical calculation of the pres- 
sure dependence of the 3rd (dotted blue line) and 5th 
(open red circles) harmonics in argon integrated over the 
full radial distribution. To be compared with the Fig. 3 
of [12]. The spectrum calculated at 50 mbar is shown in 
the inset. 



sures are 6 and 218 nJ, in good agreement with the ex- 
perimental values of 4 and 140 nJ, where losses due to 
the setup lead to a slight underestimation of the output 
energies [12]. 

If the HOKE are not considered in the model, the ra- 
tio of the FH to the TH at a pressure of 50 mbar drops 
to 0.017, and the FH and TH energies are respectively 
1.7 and 584 nJ: These values are inconsistent with the ex- 
perimental results of Kosma et al. Furthermore, contrary 
to the experimental observations [18], the FH would ex- 
hibit strong maxima at 160 and 250 mbar. These dis- 
crepancies show that the HOKE arc necessary to repro- 
duce the experimental results [12,16], further validating 
their measured values [1]. Note that the ratio of 0.017 
strongly depends on the propagation distance, so that it 
cannot be directly compared to that of 10,000 predicted 
by Kolesik et al. for the "classical" model over an un- 
specified propagation distance. For a propagation length 
of 220 /im, 80 times shorter than in our work but consis- 
tent with neglecting the phase matching, our calculation 
indeed predicts a ratio of 10,000. 

In conclusion, as recently suggested in [9], we have 
compared the recent experimental measurements of the 
TH and FH yields in argon [12] with both analytical 
and numerical simulations. These results agree quanti- 
tatively with the measured high-order Kerr indices [1]. 
This conclusion is supported by the following findings. 
First, the harmonic yield reported in argon by Kosma et 
al. at the pressure that optimized the 5th harmonic leads 
to a ratio of about 0.1 between the fifth and the third 
harmonics. This ratio implies a ratio of the Kerr indices 
consistent with our measurement of the HOKE indices 
within their uncertainty range [1] . Second, the analytical 
model based on our HOKE indices reproduces the pres- 
sure maximizing the TH, as well as the first pressure 
maximum of the 5th harmonic yield. Third, a full nu- 
merical propagation model accounting for the dispersion 
and nonlinear effects such as ionization and higher-order 



Kerr effects quantitatively reproduces the ratio of the 
harmonic yields observed in the experiment, as well as 
the pressure dependence of both the 3rd and 5th harmon- 
ics. It even reproduces the absolute harmonics intensity 
within a fairly good accuracy. 
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